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Abstract 

We propose a model for chiral polymerisation and investigate its symmetric and 
asymmetric solutions. The model has a source species which decays into left- and right- 
handed types of monomer, each of which can polymerise to form homochiral chains; 
these chains are susceptible to 'poisoning' by the opposite handed monomer. Homochiral 
polymers are assumed to influence the proportion of each type of monomer formed from 
the precursor. We show that for certain parameter values a positive feedback mechanism 
makes the symmetric steady-state solution unstable. 

The kinetics of polymer formation are then analysed in the case where the system 
starts from zero concentrations of monomer and chains. We show that following a 
long induction time, extremely large concentrations of polymers are formed for a short 
time, during this time an asymmetry introduced into the system by a random external 
perturbation may be massively amplified. The system then approaches one of the steady- 
state solutions described above. 

Keywords autocatalysis, bifurcations, chiral polymerisation, cross-catalysis, growth kinetics, 
symmetry-breaking. 



1 Introduction 

Studies of the origins of life raise many associated fundamental questions. Among these, 
one is concerned with the origin and propagation of molecular handedness. It is well known 
that chirality is a signature of life as we know it. Nucleic acids contain only D-sugars while 
proteins are made only from L-amino acids (although D-amino acides do occur in Nature 
and even occasionally show up in some proteins (Jung, 1992). What leads to the synthesis of 
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homochiral polymers, in which all the constituent monomers have the same handedness ? And 
what is responsible for the evolution of chiral purity, the more or less exclusive dominance of 
one macromolecular handedness over its mirror image? These are questions of great interest 
and importance and remain the subject of much discussion. 

It is known that, in general, the addition of the correct enantiomer to a growing polymer 
chain is more favourable that the wrong one (Joshi et al. 2000). Indeed Joyce et a/. (1984) 
showed that addition of the wrong-handed monomer to a growing oligonucleotide chain acts 
as a chain terminator, stopping all further reaction. For the case of proteins, there is also the 
driving force of the beneficial secondary structures, such as a-helices and /3-sheets, that may 
arise from homochiral polymers. Given these assumptions, that wrong-handed monomers 
inhibit chain growth, our paper is concerned with whether, starting from a racemic mixture 
of monomers, it is possible to produce a system of homochiral polymers of a greater or lesser 
degree of chiral purity. Starting from an achiral substrate, we shall be concerned with whether 
it is possible to produce a system of homochiral polymers of high chiral purity by analysing 
some plausible kinetic models. 

There is some discussion of related matters and experimental observations in the recent 
literature. Zubay (2000) provides a readable discussion of possible pre-biotic chemistry, while 
Colonna et al. (1994) describes a number of alternative self-reproducing systems. Sandars 
(2002) reviews the range and importance of chirality in biological systems, as well as the 
chemical processes which lead to achiral states. After summarising a possible historical order 
of events in the origin of life on Earth Sandars discusses the stages where a bifurcation to a 
chiral state may occur. He then applies existing knowledge of chirality on Earth to speculate 
on the question of extra-terrestrial life and its chirality. 

Luisi's group at ETH in Zurich has studied various polymeric systems in which left- and 
right-handed monomers aggregate together to form larger than expected concentrations of 
homochiral polymers (Blocher et al. 2001, and Hitz et a/.2001). Hitz et al. present an analysis 
of the data and postulate that the exess of homochiral polymers is due to a high-order Markov 
process rather than the feedback mechanism which we analyse in this paper. It is readily un- 
derstood that the rate coefficient governing polymer growth may depend on the handedness of 
the monomer to be attached to the chain and handedness of the monomer which currently ter- 
minates the chain; however, with a high-order Markov processes this rate coefficient may also 
depend on the handedness of the penultimate (and possibly the antepenultimate) monomers 
in the polymeric chain. 

While ohgopeptides spontaneously form homochiral sequences, Hitz & Luisi (2002) have 
shown that the presence of quartz promotes the production of a high yield of homochiral 
sequences. More recently (Hitz & Luisi, 2003) this has been quantified and the level of 
enantiomeric exess in a system has been amplified from 20% to over 70% and in some cases 
to 100%, by the presence of quartz. 

The model of symmetry-breaking in chiral polymerisation which we explore in this paper 
is based on a model suggested by Sandars (2003). Sandars gives an account of the history 
of chemical discoveries leading up to the mechanisms of enantiomeric cross-inhibition and 
autocatalysis upon which his and our models rely. Sandars integrates the resulting system of 
ordinary differential equations numerically to explore the parameter regimes where symmetry- 
breaking solutions exist. He observes a 'phase-transition' type of phenomenon where a small 
change in the fidelity of the feedback mechanism leads to a large-scale change in the steady- 
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state which the system as a whole converges to. Below a critical fidelity in the nonlinear 
feedback process the system approaches a symmetric state where equal amounts of left- and 
right-handed polymers coexist, whilst above the critical fidelity a homochiral state is ap- 
proached in which one chirality of polymer dominates to the almost complete exclusion of the 
other. 

The model studied here includes the inhibition of homochiral sequences of long chains. 
This bears some similarity with our modelling of cement hydration (Wattis & Coveney, 1997) 
in which larger clusters were susceptible to poisoning by another component. Here the poison 
is simply the monomer of opposite chirality, so each monomer can either act gent of 

growth (of polymers of the same handedness) or of inhibition (of polymers of the opposite 
handedness). This dual role leads to some subtle and interesting effects since an abundance 
of, say, right-handed polymers makes it unlikely for any left-handed polymeric sequences to 
form, and the majority of left-handed monomers produced will be consumed by inhibiting 
right-handed homochiral sequences - a 'double- whammy' effect. This form of competition is 
distinct to the models of nucleation involving competition we have analysed previously; for 
example in Wattis (1999) and in Bolton & Wattis (2004) there is only one monomer which 
assembles to form two morphologies of cluster. Competition is thus between the growth of 
one type of cluster and that of another. However, similar mechanisms are operative in both 
those examples and in the present paper, since both types of homochiral polymer sequence 
(left and right) are ultimately composed from, and hence competing for, the same source 
material. 

We have investigated the growth of RNA chains in an earlier paper, (Wattis & Coveney, 
1999) wherein we used a much more detailed and hence complicated model to assess the feasi- 
bility of long self-replicating RNA sequences forming within a realistic timescale. While that 
model of RNA polymerisation had no precursor species, it contained four types of nucleotide 
monomer (A, C, G and T), both autocatalytic and cross-catalytic polymerisation mechanisms 
and an important hydrolysis step that recycled growing RNA sequences. On the basis of plau- 
sible assumptions about the prebiotic soup, and by invoking a number of approximations and 
coarse-graining procedures, we were able to show that self-replicating RNA sequences are 
amplified in such mixtures, while less capable replicators are driven to extinction. 

By comparison, in the present paper we are able to construct a model which is more di- 
rectly tractable using standard methods of mathematical analysis. The paper is structured as 
follows. In section|2l we specify the basic kinetic model, while in section|21we study its steady- 
state solutions. Section|3]considers the time-dependent achiral solutions and investigates their 
kinetic stability. Section El then considers the case of perfect chiral symmetry-breaking. It is 
followed by a discussion (section IHl) and conclusions from our work are presented in section [7| 

2 The kinetic model 
2.1 Microscopic modelling 

We aim to investigate spontaneous symmetry breaking in a system which allows both right- 
and left-handed chiral polymers to form. We assume there is some achiral source S which 
spontaneously transforms into right- and left-handed monomers at a slow rate e; we also 
assume that the presence of longer chiral polymers (denoted L or R) accelerates the formation 
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of monomers of the same chirality. Thus we shall be concerned with the following set of coupled 
chemical reactions: 

HHf) Hi-f) 

S^L,, S + Q ^ L^ + Q, S + P^L, + P, 
S^Ri, S + P^Ri + P, S + Q ^ Ri + Q, 

where Li, Ri are the left and right monomer species respectively while Q (P) represents some 
measure of the total concentrations of left-handed (right-handed) homochiral polymers in the 
system. The precise forms of these rate processes will be specified later (see equation (jH)). 
The rate coefficients of the reactions are e, k. We shall assume that e <C fc, and that k depends 
on the length of the polymer L, R. The parameter / in these rate coefficients is the fidelity of 
the feedback mechanism; typically this will not be perfect, that is, / < 1 is likely in general. 

This model has some similarities with models we have studied in other areas of investi- 
gation; for example, the kinetics of micelle-formation in ethyl caprylate (Coveney & Wattis, 
1996). In this system the breakdown of caprylate ester into monomer occurs spontaneously 
at some slow rate, but is massively accelerated by the presence of micelles, which have a 
catalytic role in the breakdown of the source species. In the present paper this mechanism 
is more complex since there are two monomers, and there is an additional fidelity parameter 
since long left- or right-handed sequences can promote the formation of the oppositely-handed 
(right/left) monomer as well as its own. In the caprylate system, unusual kinetic behaviour is 
observed when the system is initiated without any product present. The system is then effec- 
tively in a metastable state, and very little appears to happen for a long time. Following this 
induction time, the kinetics then proceed fairly rapidly. We expect to see similar behaviour 
in the current model in the case e <^1. 

Sandars' model (2003) differs from ours in that his imposes a maximum polymer length, 
typically set at five, and only polymers of this maximum length act catalytically in the 
breakdown of the source into monomers. We allow polymers to grow to arbitrary lengths, 
and chiral polymers of all lengths have some degree of efficacy in the autocatalytic feedback 
mechanism by which the source species decays to form chiral monomers. 

The monomers will be allowed to combine to form chirally pure polymers, denoted by L„ 
and Rn according to 

Ln + Li s> Ln+l, Rn + -Rl ^ Rn+1- (2) 

We assume that the monomer of opposite handedness may attach to a growing polymer and 
so inhibit its growth; such inhibited sequences will be denoted by RL^ for a polymer L„ which 
has been terminated by an Ri monomer, and LR^ for the corresponding -R„Li polymer. We 
denote the rate of such reactions by ax- Thus we have the two rate processes 

Ln + Rl RLn, Rn + -^1 LRn- (3) 

However, the sequences RLni LRn are treated as inert products which have no influence on 
the other rate processes; therefore their concentrations can be ignored in the mathematical 
modelling of the chemical reactions. Following Joyce et al. (1984) we assume no further 
growth of these products can occur. The system studied by Joshi et al. (2000) corresponds 
to < X < 1; however, we shall consider the full range of possible x > 0- 
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2.2 Macroscopic model 

To close the system of rate equations ensuing from the scheme (^-(jBI) we shall assume that the 
source chemical, S, is added to the system at some constant rate 5*0. We define macroscopic 
quantities for the total concentration of all homochiral sequences of each type of polymer by 
L and R, and the mass of monomers in each set of homochiral sequences by Q, P in schema 
©by 

oo oo oo oo 

L = 5]L„, R=J2R^, Q = Y^nL^, P=Y^nRn. (4) 

n=2 n=2 n=2 n=2 

Applying the law of mass action to (H}-© and using the definitions of (jU, we obtain the 
kinetic equations 



dS 

Itt 
dLi 

IT 

"df 
dRi 

"df 

dRn 



So - 2eS - kS{P + Q) (5) 

eS + \kS[{l + f)Q + (1 - f)P] - aLi(2Li + L + xR) (6) 

aLi(L„_i - Ln) - axLnRi {n > 2) (7) 

eS + lkS[{l + f)P + (1 - f)Q] - ai?i(2i?i + R + xL) (8) 

aRi {Rn-i -Rn)- axRuLi (n > 2) . (9) 



This constitutes a quantitative description of the polymerising system. Equations ((Tj) and © 
hold for all integers n > 2, that is for all polymer lengths, from dimers (n = 2) to infinitely 
long polymers. A diagrammatic summary of the rate processes involved is shown in Figure ^ 
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Figure 1: Diagrammatic representation of the homochiral polymerisation scheme under study. 

A considerable advantage of the present model is that it is possible to reduce the complexity 
of the corresponding infinite set of rate equations ©^Q? resulting in a closed system of only 
seven equations which contains the dynamics of the full system. This is based on the quantities 
S, Li,Ri, L, R, P, Q which evolve according to 

— = So-26S -kS{P + Q) (10) 
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~df 
dL 

It 
dQ 

dt 
dRi 

"df 
dR 

dP 

Itt 



eS + lkS[{l + f)Q + (1 - f)P] - aLi{2Li + L + xR) (11) 

aLl - axLRi (12) 

2aLl + aLLi - axQRi (13) 

eS + lkS[{l + f)P + (1 - f)Q] - aRi{2Ri + R + xL) (14) 

aRl - axRLi (15) 

2aRl + aRRi - axPLi. (16) 



The rate processes described in terms of equations (fTnjl - lfTHIl are depicted in Figure El 
This reduction from an infinite system of coupled ordinary differential equations to a system 
of just seven ordinary differential equations is remarkable in that it is an exact simplification 
- no approximations have been made. Such an exact reduction would not be possible if the 
growth rate coefficients were dependent on polymer length. If appropriate at all it would 
then have to rely on approximations, and an accurate approximation may require not just 
the zeroth and first moments of the distributions i?„, but on higher moments as well. 

e + lk{l + f)Q 




e + \k{l + f) 



Figure 2: Illustration of our homochiral polymerisation scheme including the nonlinear feed- 
back mechanisms. The quantities LR^, RLn are not included in the mathematical model 
since they play no role in the reaction scheme. The quantities L, R refer to the total numer 
of polymers of each handedness, and P, Q to the total mass of material in polymeric form. It 
is these latter quantities (P, Q) which determine the effectiveness of the catalytic breakdown 
of the source species into monomers. 



2.3 Transforming the system of kinetic equations 

Before we analyse the steady-state kinetics, it is useful to recast the system of seven ordinary 
differential equations in an alternative form. Instead of describing the concentration of each 
species separately, we assign variables for the total concentrations of source material (S), 
monomers (/i), sequences (A^) and sequence mass (M), 

fi = Li + Ri, N = L + R, M = P + Q, (17) 



6 



and a set of variables S, 6, t] for the proportions of right-handed molecules in each of /i, A^, M 
respectively 

R,-L, R-L P-Q 

These quantities relate to the chiral purity of the system: 5 is the chiral purity of the monomers 
(yu), 9 and 7] describe the chiral purity of the homochiral polymer chains, 9 being a number- 
weighted measure (corresponding to A^) and rj a mass-weighted measure (corresponding to 
M). 

We transform the kinetic equations ()10 p - (fT^ into the variables relating to the concentra- 
tions of polymers given by (fT7j) - (fTHj) and hence obtain the seven coupled equations 

-- So-2eS -kSM, (19) 



dt 

^ = 2eS + kSM-afj,^il + 5^)-^aixN{l + S9)-^axi^N{l 



(20) 



dN 

'dt 
dM 

~~dt 



\a^l\l + 5^)-\ax^lN{l-59), (21) 
a/i^(l + 5^) + ia/iiV(l + 59) - \axiiM{l - 5ri), (22) 



f = ^i2S-v-v6') + ^{6 + 9-r^-69r^)+',axf.S{l-v% 

(23) 

f = ^{26 - 9 - 96') + laxf^6{l-9'), (24) 

d5 2eS6 kSM6 kfSMr] , ,^ ^.n.o . 

— = + — --la{l-5^){2i^6 + N9-xN9). 

dt n n jj, ^ 

(25) 

The advantage of (fT!Hl - (j^ over (fTn|l - (pi)|l lies in the considerable reduction in the amount of 
algebra required to derive solutions; for example it is easy to see that S = 9 = rj = satisfies 
leaving a system of four equations (fT^ - (j^ for the four unknowns S, /i, A^, M. 
The equations (fTUjl - lj^ will be used throughout the rest of the paper. In later sections 
we shall consider the equations for the symmetric growth of homochiral polymer sequences, 
(fnijl - lj^ . separately from the equations describing the chiral purity of the system, (f^ - pHjl . 



3 Steady-state behaviour 

For the analysis of the steady-state solutions all the right-hand sides of equations (fTUIl - lfTHI) . or 
equations (fT^ -(P3 |) are set to zero. For the later stages of the calculation of the steady-states 
we ignore the small parameter e (that is, we set e = 0), since this will allow explicit analytical 
formulae to be derived and this will not greatly influence the steady-states. The reason for 
this is that once there are appreciable numbers of polymers of either handedness present in 
the system, the catalytic breakdown of source will dominate the production of monomers of 
both handednesses. The 0{e) terms in equations p()|) - ()lfij) will be reinstated later where the 
kinetics of the system starting from zero initial data are investigated. 
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We solve the equations (fTUj). (fT^ and (fT3j) to express the solution in terms of Li, Ri 



. Ll Rl Rl{R, + 2xU) Ll{U + 2xRi) 

xRi xLi x^Lf x^Ri 

2ex^L\Rl + k{L\ + R\ + 2xL\Ri + 2xLiR\)' ^ ' 

The microscopic steady-state solution is then found to be 

In terms of the variables N , M, 6, 77, the steady-state solution is then 

_ Ml + 3^^) , . _ /x[l + 10^^ + 5^^ + 2x(l-^^)(l + 3^^)] 



together with 



2£x'(l- 52)2 + + 1052 + + 2x(l-52)(l + 352)]- ^ ^ 



The average chain length is given by M/N, which can be written as 

M 1 + 10^2^5^4 

iv x(i-<5')(i + 352)- ) 

The most important thing to notice about the second term here is that asymmetric (chiral, 
i.e. S ^ 0) solutions permit longer chains to be produced. For the symmetric solution (5 = 0) 
we have M/N = 2 + 1/x', whereas if S approaches ±1 then arbitrarily large chains can be 
produced. We note that large inhibition rates (x) reduce the expected sequence length (as 
one might expect); however, we shall show later on that this effect is compensated for since 
larger values for x make it easier for the system to adopt a chiral state. 

It is more natural to express solutions in terms of fi, 6 than Li, Ri. The chiral purity of 
the polymers can be defined in two ways, one by the mass-weighted purity (77) and the other 
in terms of the number- weighted purity (6). At steady-state we have 

^ . (1±^\ , = . / 5 + 10^2^^4^2x(l-^^)(3 + ^^) \ 

1^1 + 352^' ' 1^1 + 1052 + 5(^4 + 2x(l- 52) (1 + 352) y ■ ^ ^ 

Thus if the chiral purity of the monomer, 6, departs from zero, then the chiral purity of the 
homochiral sequences {9 and 77) also do, and by greater amounts, as illustrated in Figure El 
This can be seen from the linearisations for small 6 which give 6 = 36 and rj = [5 + 2%) 5/(1 + 
2x). For large x the curve for rj approaches that for 6. 

When substituted into equation (j2(Jj) the above equations yield 

^2 = ^xSo{l-P) 

^ 3ax(l - 5^) + a(l + 652 + 54) • ^ ^ 
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Figure 3: Graph of steady-state values for the chiral purities {9, if) of the homochiral sequences 
against the chiral purity of the monomer, 5. These are given by equation (jH^ . and show 
the accentuated effect of symmetry-breaking in the homochiral sequences compared to the 
monomers. The solid line corresponds to 9; the dashed lines refer to rj. The steepest curve 
relates to x = 0, the next to x = 1 and that for x = 2 is the closest to the 9 curve. 



Finally we have to determine the chiral purity of the monomers, 5. This is the most compli- 
cated part of the calculation so, for the sake of clarity, we now set £ = 0; 5 is then given by 
either 5 = 0, the symmetric (achiral) steady-state solution or, in terms of the fidelity, by 

^ 1^ 5 + 10(52+54 + 2x(l-52)(3+52) j l^i + 652 + 54 + 3x(l-5^)y ■ 

We shall discuss this chiral solution in further detail later on (Section 13. 2j) . One question we 
aim to address in the remainder of this section is how small / could be, and the system still 
exhibit a bifurcation to a non-symmetric state. 



3.1 The symmetric, achiral, steady-state solution 

It is not immediately obvious that all the solutions given above in equations (P7j) - (jM|) for 
the steady-state exist for all parameter values, or whether any particular solution is unique. 
Hence we start by considering the symmetric solution, where 5 = ^ = r/ = 0. This solution 
exists for all parameter values and in this case we have 



/ 2x^0 ^ xVqX^o(1 + 3x) 

^"V«(1 + 3X)' V2k{l + 2x) ' ^ ^ 



250 ^ M=il±^ = ii;M^, (36) 
X V«X(1 + 3X) Xv/«X(1 + 3X) 

with Li = Ri = \yi,L = R=\N,P = Q = \M. 

We now consider the mathematical stability of this solution, that is we aim to answer the 
question: 'If the system is close to the steady-state solution, will it be attracted closer to it. 
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Figure 4: Graph of the chiral purity S against x ^-^id / for the steady-state values given by 
equation for fidehty in the range < / < 1 and with < x < 2. 



or diverge further away from it ?' To answer this question, we hnearise around the solution 
5 = 6 = 7] = 0. Note that the formulae for fi, N, M, S all have 0{6'^) correction terms, and no 
0{6) terms. Thus these will be treated as constants. We only need to analyse the evolution of 
{6, 6, rj) over time, starting from small perturbations away from (0,0,0). Linearising equations 
|)- (f^ we obtain 



drj 
lit 

de 

d6 



M 



2M 



{25 -e) + iax/i5 
kfSMr] 



2N 
kSM6 



+ 



- \a{2^i5 + Ne-xNe). 



(37) 
(38) 
(39) 



Note that we have ignored the 0{e) term. Inserting the expressions for /x, and M from 
(123)-® we obtain 



d 
It 





( 


(:! 











/(l+3x) 




(40) 



which has the characteristic polynomial 

x'X' + {2x' + 5x + l)X' + 







X 



7x + 5 



2/(1 + 3x)- 
l + 2x 



A + 



+ 4 + 2x - 



/(5 + 6x)(l + 3x) ' 
(1 + 2x) 



(41) 



The Routh-Hourwitz criteria state that all solutions of X^ + AX^ + BX + C = satisfy 3fJ(A) < 
if and only if A > 0, C > and AB > C (for further details, see Murray, 1989). If all values 
of 3?A are negative then any solution of (pn|) will have all the quantities rj, 9 and 5 decaying to 
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zero as time increases. In our example, clearly A > whatever values x ^i-nd / take, but the 
other two conditions are less clear. The condition C > gives / > fdx) with fdx) given by 



2(2 + x-)(l + 2x) ,,„, 
^=<''' = (5 + 6x)(l + 3x)- 



This value for fc agrees with / as given by ()34|) in the case 6 = 0. The condition AB > C 
implies an instability when / > fc2{x)'i however, fc2{x) liss in the region / > 1 for all x so 
this instability can be ignored, since only / < 1 is physically realisable. So the symmetric 
solution is stable for / < fc, and unstable for / > fc. The value of fc lies between 2/9 and 
4/5 and depends on x- at x = 0; the achiral solution is unstable for / > 4/5, whereas at large 
X, the instability of the achiral solution occurs for / > 2/9. 



3.2 The chiral (asymmetric) steady-state solution 

Once / > fdx) we have the existence of two chiral steady-states as well as the achiral steady- 
state. Furthermore, since the Routh-Hourwitz criteria are 'necessary and sufficient' for the 
existence of a solution of (PT|) with 3?(A) < 0, once f > fc the achiral solution becomes 
unstable, so will not be observed in any physical system. Thus once the chiral solutions exist, 
they become the 'preferred' steady-state solutions since they are mathematical stable. Even 
if a system were artificially put into the symmetric (achiral) state, any small perturbation 
would cause some chiral imbalance and the natural kinetics of the model would then carry the 
whole system to one of the two steady chiral states in which one handedness of homochiral 
polymers dominates the other. 

At / = /c a bifurcation occurs, and two mirror-image chiral solutions appear. We now 
briefiy examine the neighbourhood of this point in more detail. These chiral steady-states 
are governed by equation (j2D); from this equation we note 6 = implies / = fc- Assuming 
small 6 and expanding we find f = fc + 136"^ with /? > so the bifurcation is supercritical, and 
the new solutions exist in f > fc where the achiral solution is unstable. Thus, as expected, 
the bifurcation occurs at exactly the same position as that at which the achiral steady-state 
becomes unstable. This is a standard result in bifurcation theory; see Guckenheimer & Holmes 
(1983) or Berge et al. (1984) for more details. Although this is a curve in (/, x) space, we shall 
consider / to be the primary bifurcation parameter, with the bifurcation point fc depending 
on X- 

In figure m we show the chiral steady-state solutions given by (jMj) parametrically, plotting 
6 against / and x- The bifurcation occurs at fc = 4/5 if x = 0, reducing to fc = 2/9 in the 
limit of large x- Thus we see the beneficial effect of sequence inhibition if one wants a system 
which undergoes a symmetry-breaking bifurcation at small values of the fidelity parameter /. 
Figure El shows the chiral purities of the polymers, firstly weighted by number {9) and then 
by mass {rj). Both show that the chiral purity of the polymers is much greater than that 
of the monomer (compare Figure El with Figure E)). Note also that for the range of fidelities 
2/9 < / < 4/5 the bifurcation to asymmetric (chiral) solutions can occur by increasing x- 

Another natural question to pose here is whether only small and moderate values of 6 can 
be accessed; we put v = 1 — 6 with v <^lio determine under what conditions 6 can approach 
unity. We find from equation that f = 1 — xi^- In this case we also expect 9 and rj to 
be near unity. In fact the chiral purity of the polymeric sequences is much enhanced over the 
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Figure 5: Graph of the chiral purities 9 (on the left) and rj (right) against x ^i-nd / for the 
steady-state values given by equation (jSlj). This is to be compared with Figure El illustrat- 
ing that the polymer concentration and mass show more extreme chiral purities than the 
monomers (described by 5). 



chiral purity of monomers, since a two-term expansion of equation ()32|) gives 

e^l- \v\ ~ 1 - \xv\ (43) 

for z/ -C 1 (in the special case x = we have ~ 1 — ^z^^). The concentrations scale with v 
according to 

V 2a V (^X^ V O'X 



^^-m ^^m 

together with S = 0{y^''^). 

In summary we find that chiral solutions do indeed exist provided f > fc and this is easier 
to satisfy at larger inhibition rates (x)- Thus the presence of stronger cross inhibition aids 
the manifestation of chiral steady-states. Also, the higher the fidelity (/), the greater the 
dominance of one chirality over the other. 

We now turn to an analysis of the kinetics of sequence growth, in order to determine at 
what stage of the reaction the system is likely to manifest a chiral state. 



4 Kinetics and stability of achiral solutions 

In this section we reintroduce the 0{e) term omitted in the analysis following (|32p and solve 
the kinetic problem for the achiral solution, in the limit £ ^ 1. We assume that initially there 
are no monomers (yU = 0), no polymers (M = N = 0) and no precursor (5 = 0), though this 
source material is added continuously to the system starting at time t = 0. If e = then since 
there are no polymers, the precursor cannot be broken down and no polymers will ever form, 
so we need the 0{e) term to produce some monomer. In section we analyse the stability 
of the solution, that is, whether small, random perturbations to such a solution grow or are 
damped out, as the total number of polymers and monomers grow from zero concentrations. 
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4.1 Growth of the achiral solution 



We assume a set of initial conditions in which rj = 6 = = 6, and then rj = 6 = = 5 for all 
subsequent times. This reduces the system (fTI?|) - (f^ to a system of three ordinary differential 
equations which we solve asymptotically in the limit £ -C 1. Our aim is to develop matched 
asymptotic expansions for the solution of 



dS_ 

H 

dji 

It 

dN 

~dt 

dM 

through a series of timescales. 



So - 2eS - kSM, (46) 

2eS + kSM-afj.'^-lafxN{l + x), (47) 

ia/i^ - ^afixN, (48) 

a/i^ + ^aNfi - \axM^, (49) 



4.1.1 Timescale I: t = 0{6-^/^) 

For this timescale the appropriate scalings are t = e^^^^ti together with 

S = e^^/^Su /i = 5=^/^1, N = eNi, M = eMi, (50) 
thus this is a long induction time over which the leading order equations are 

S[ = So, /i'l = 2Si + kSiMi, N[ = ^afil, M[ = afij, (51) 
where prime denotes d/dti. This system has the solution 



Si = Soh, N, = lM,, Mi = -|i---. (52) 
^ kSoti k 



where fii is given by the solution of 



,1 



fi'l - ^ - Sokahul = 0, (53) 
^1 

which unfortunately is not explicitly available. However, we can see that the timescale ends 
abruptly with f^i,Ni, Mi all diverging as ti tic, according to 

6 6 , , 

akSotUtic-tiY' ak^SltUtic-tir ^ ' 

for some constant tic. These relationships help us determine the scalings relevant in the next 
timescale. In this timescale we have seen the accumulation of source material, but this is 
only slowly converted into monomers and chains, so both of these grow very slowly, causing 
a big build up of source material until, at the end of this timescale, we see the concentration 
of chains increase to the level where the catalytic mechanism becomes active and accelerates 
the formation of monomers and chains. 
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4.1.2 Timescale II: t = e'^'Hi^ + 0{e^/^) 

In this timescale all quantities are large and evolve quickly. To be specific we have 

S = e-^''>S2, ^ = £-1/5^2, N = e-^'''N2, M = ^-^/^Ma, (55) 

together with t = e~^/^tic + e^^^t2. Using primes to denote d/dt2, the leading order equations 
are 

S'^ = -kS2M2, fi'2 = kS2M2 - afil - ia(l + x)l^2N2, 

N2 = \aix\ - iax/i2A^2, M'^ = a^l + \a^2N2 - |ax/^2M2, 

(56) 

As the chains are present in large enough quantities for the catalytic mechanism to be active, 
and since there is a large amount of source material present at the start of this timescale, 
this source material is rapidly converted into monomers and chains so that /i2, ^'"2 and M2 
all increase at the expense of the source species, whose concentration now monotonically 
reduces, so that the only significant simplification in the equations is in the equation for 5*2. 

No explicit solution which matches back into Timescale I is available; however, the form 
of the large-time solution in Timescale II can be determined. Consider new timescale given 
by — ^ = ^, hence ^2 = f / "7^)5 then the system can be written as 



/J2 dt2 dr'' ^ a J ii2(f)' 

dH2 

dr 
dN2 



Q - 2/i2 - (1 + x)N2 (57) 
/i2 - XN2 (58) 



dr 

dM2 , , 

= 2fX2 + N2-xM2, (59) 

which is a linear system, together with the equation 

— log52 = , 60 

dr a/i2 

from which Q = 2kS2M2/afi2 is obtained. The solution of equations ()57|) - (j3^ is given 
by a combination of a complimentary function which solves the system with Q{t) =0 and a 
particular solution which satisfies the Q{t) input term in equation ()57|1 . For general parameter 
values, the forms of these solutions cannot be explicitly determined, however we shall study 
the two particular limits of large and small x in which approximations can be obtained. These 
correspond respectively to the cases where the solution is dominated by the particular solution 
and the complimentary function at large times. 



4.1.3 The solution in timescale II for x ^ 1 

If we assume that the input function Q has the form Qoe"^'^, then the particular solution has 
the form 

/ /^2(r) \ ( fi2 \ 

N2{r) \ = \ N2 e-^\ (61) 
\M2{t) J \mJ 
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Figure 6: Plot of a numerical solution of 5", /x, N and M against time for the case e = 10~^, 
a = 1, k = 0.8, 5*0 = 1, X = 3.333. The dotted curve corresponds to S(t), the dsashed curve 
to fi(t) and the upper solid curve to M{t) and the lower solid curve to N{t). 



In the calculation of Q, we then have 5*2 



-Ar 



too, and this assumption imphes A 



2kM2/afi2, and Q ~ AS'2e~'^'^. In the calculation of the prefactors fi2,N2, M2 we then have to 
solve a cubic, and there are constraints that all the prefactors must be positive. 
Inserting these into the equations (j37|) - (j^ we obtain 



1^2 



/i2(2x-2A + l) 



X 



A 



and the cubic equation for A is 



X'-2x\'+ x' + 



Ak 



A- — l + 2x =0. 

a 



(62) 



(63) 



For asymptotically large X) this has two roots near \ = x-, and one near A = 0. The larger 
roots are at A = x =t ^J^mJqtx-, one of which violates the condition A < x and the other leads 
to a solution which rapidly decays in time. The physically relevant solution corresponds to 
A ~ Ak/aXi hence the solution 



(52,/i2,iV2,M2) ~ A 



Ak 



X,l,2 e 



-Akr/ax 



(64) 



Since this result has been derived on the basis of large Xi we see a slow exponential decay in 
r. The complimentary function also decays exponentially in r, but with exponents Ai = 3, 
A2 = X ~ 1 and A3 = Xi hence the complimentary function decays much more rapidly than 
the particular solution. 

In terms of the original timescale ^2 the solution (jMj) leads to 



(52,/i2,Ar2,M2) 



'3ax^ X 1 r 
%kH2 ' 2A;t2 ' 2H2 ' kt2 



(65) 



This agrees with the numerical observed results, which suggest that for x ^ 1 we observe that 
all quantities decay with 1/^2- Figure shows a numerical solution of the system produced by 



15 



Matlab 6.5.0 Release 13; the right-hand graph shows in detail the decay of the concentrations 
at the end of the second timescale. When the concentrations S,n,N,M become 0{1) then 
other terms become relevant in the kinetic equations, and a further timescale is required to 
describe the final approach to equilibrium. 

4.1.4 Timescale III for larger x- t = ^~^^^hc + 

As 5*2, N2 and M2 all decay like l/t2 at the end of the previous timescale, the new first 
term to become significant is the 5*0 input into the S equation. This becomes significant when 
^2 = e"^''^, thus the third timescale is t = £:~^/^tic + ^3. In this timescale all of 5, /i, M and 
are 0{1). Thus the leading order equations are 



dS_ 

'dt 
dfi 

'dt 

dN 

"dt 

dM 

"df 



So - kSM, 

kSM - a/i^ - |a/iA^(l + x), 



(66) 
(67) 
(68) 
(69) 



and over this timescale the system approaches its steady-state. Since at leading order the 
eS terms are neglected, the leading order steady-state solution approached is precisely that 
described in Section |21 

Overall, we see that there is a long induction time, of 0{e~^^^), followed by a some rapid 
kinetics during which the system explores states a long way from its steady-state, and then 
over a relatively fast (C(l)) timescale the system approaches its steady-state. We now derive 
the kinetics of the achiral solution for smaller x before going to address the stability of the 
growing achiral solution. 



4.1.5 The solution in timescale II for x ^ 1 

Numerical simulations suggest that at the end of TII, 5*2 and ^2 decay exponentially (in ^2) to 
zero, with M2 and A^2 tending to constant values; eventually reaching their C(l) steady-state 
values over subsequent and much longer timescales (see later subsections and Figure [7| for 
example) . 

The linear system (jFTjl - ffHUI) is solved by the sum of a particular solution and a compli- 
mentary function. For smaller x the large-time solution is dominated by the complimentary 
solution component of the solution (i.e. that with Q = 0) which is given by a combination 
of exponentials. In the case x the eigenvalues are Ai = — x, A2,3 = — 1 ± iy/2x, thus 

applying the special result for repeated roots, at leading order, we have the complimentary 
function 




/o 





/ 



C2 




(70) 
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Figure 7: Plot of numerical solution of S{t) (the dotted curve), (the dashed curve), N(t) 
(the lower solid curve) and M{t) (the upper solid curve). All parameters are as in Figure [7| 
except X = 0.333. 

The solution in this case corresponds to the scenario in which /i2 — > as r ^ Tc < oo. Let us 
assume 

/i2 ~7l2(T-c-r), N2^N2, M2^M2, as r^r„ (71) 
then ^2 = f / j^^:;^ implies Tc — r = ^-"4^2^2/2 w^qx r ^ corresponds to ^2 00. We thus 
have jU2 ~ 7l2e~"'^2*2/2_ ^^ie solution of equation (plj) yields 

^2 ~ S2{rc - r)2'=^2/"^2 = 526"'"^^*^ (72) 

If ^a]l2 > kM2, then this solution predicts that fi2 decays faster than 5*2; however, towards 
the end of the second timescale the dominant terms in the equation for /X2 are /ig = kS2M2 — 
^afi2N2- This equation implies that /i2 cannot decay faster than S2, and thus we expect 
|a7l2 < kM2- If equality holds, then both /x and S reach 0{e^^^) at the same time. In 
summary, towards the end of this timescale, we find fj,2 and 5*2 decaying exponentially in ^2 

with 

^2 ~ S2e-''^''\ ^I2 7i2e-"^2*^/^ (73) 

together with N2 — > N2 and M2 — > M2. This agrees with our observations of the numerical 
solution discussed in the opening paragraph of this subsection. The system then passes 
straight into Timescale III, which becomes relevant when kSM ~ So- This occurs when 
S = C(ei/5) ^hich happens when t2 ~ t2c = (2/(5fcM2)) log(l/e). 

However, if kM2 > \(ip2 then 5*2 decays at a faster rate than yU2 and there is an further 
timescale between this and timescale III (below) over which S saturates while /i decreases 
further. This timescale is given by 

S = e'^'^S, ^i = e^fi, N = e-^'^N, M = e'^/^M, (74) 

for some —1/5 < 7 < 1/5 and t = e'^^^^tic + 6^^^^\og{l/e)t2c + s^^^t. The governing equations 
are then 

S' = Sq- kSM, p.' = -\a0, N' = 0, M' = 0, (75) 

17 



where prime denotes d/dt. Over this timescale, M and N remain constant, S equihbrates to 
So/kM and /I continues to decrease exponentially, with rate —^aN. When /i reaches 0{e^^^), 
timescale III is entered; this occurs when t = 0{log{l/s)), so the relationship between ^3 and 
t in Timescale III given below remains valid following a redefinition of t2c to incorporate this 
extra shift in time. 

4.1.6 Timescale III for smaller x- t = £~'^^%c + e^l^ log(l/£)t2c + 0{s}l^) 
In the new timescale we have 

5 = £^/%, /x = £^/V3, iV = £-^/5iV3, M = £-i/5M3, (76) 

together with t = 1^8'^/^ + tac^^^^ log(l/e) + Using prime to denote time derivative 

with respect to the new time variable, ^3, the leading order equations are 

//'3 = kSsM^-la{l + x)l^3N^, = 0. ^" ^ 

Over this rapid timescale, M and N do not change from their values at the end of the TII 
whilst // and S equilibrate to the values 

'S'n 25*0 , s 

TTT' *3 OO- (78) 



kM^' a(l + x)iV3' 

The evolution of M and N occurs over a longer timescale with ji and S constrained to their 
respective local equilibrium values. 

4.1.7 Timescale IV for smaller t = e'^^Hic + C(£^/^) 

Since all concentrations approached constants at the end of Till, their magnitudes remain 
unchanged for the fourth timescale, only the scaling for t changes, we now write 

S^e^/^S^, i_i^e^/^i_i^, N^e-^/^N^, M^e'^'^M^, (79) 

with t — ticS'^/^ + t2c£^^^ log(l/e) + e~^/^t4, and so obtain the equations 

= So-kS^M^ = kS4M4-la{l + x)fi4N4 s 
Ni = -^axfi^N^ M'^ = ia/^4(iV4-xM4). ^ ' 

These imply kS^Mi — So and IJL4N4 — 2So/a{l + x), thus we have the solution 

1 + X ' ^^-ax(t4.-t4)' ^'-kM,^ ^^'^ 
M4 = CiU, - U) - ^o(^^c-t4)log(t4.-t4) ^ ^g2) 

for some constant C. This timescale ends with yLi4 and 5*4 increasing hence becoming larger 
in the next timescale, and M4, iV4 decaying hence being smaller in the next timescale. 
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4.1.8 Timescale V for smaller x- t = ^^^^^ihc + he) + 0{l) 

The final timescale is given by all of 5*, fi, N and M being 0{1) and varying on an 0{1) 
timescale, thus corresponds to the approach to the global equilibrium solution given by (jH^jl - 



4.1.9 Summary 

For all values of the parameters, the solution starts with a long induction period during which 
the concentration of the precursor species 5* becomes large. There follows a short period of very 
rapid kinetics where all concentrations become large, and then the precursor and monomer 
concentrations decay. For larger values of x the concentration of polymers also decays and 
the steady-state is reached relatively rapidly. For smaller x the monomer concentration and 
that of the source species become very small and the polymer concentrations (mass-weighted 
and number weighted) both remain high and slowly evolve to their steady-state values over 
a longer timescale, which is of similar length to the induction timescale. These two distinct 
behaviours are illustrated in Figures El and 



4.2 Stability of the achiral growing solution 

Having determined the form of the kinetic behaviour for the achiral solution [r] = 6 = 6 = 0), 
we now consider the linear stability of this solution. Assuming there is some small random 
perturbation during the evolution, we ask whether a perturbation grows or decays as time 
progresses. We shall use the already determined solution for S{t), fi{t) N{t) and M{t) (from 
section ITTj) . and assume that the perturbation does not make any alteration to these total 
concentrations of source, monomer and polymer at leading order. This assumption was seen to 
be valid in the case of the steady-state solution, and since the kinetic equations are symmetric 
under the transformation (5, 9, rf) (—5, —6*, —rf), we expect modifications to /i, M, S also 
to be of second order in S, 9, rj. From equations we have 
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4.2.1 Stability of achiral solution in Timescale I 

We focus our attention on the linear stability of 6(t), 9{t) and rj{t) as given by In 
the first timescale, using the scalings (j50|) . we find the simplified linear stability problem 
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2kSiMi 
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) 



It is clear that this matrix has a simpler structure than the general case, which has only one 
zero entry, as given in equation ()4()j) . The eigenvalues of the matrix above are all negative if 
/ < 1/2, but if / > 1/2 then one is positive, indicating that a perturbation away from 9,rj,5 = 
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would increase in size as time progresses. The temporal evolution of such perturbations is 
non-trivial since Mi, Ni, fii and 5*1 are all time-dependent. 



4.2.2 Stability of achiral solution in Timescale II, larger x 





From equation ()83|) for the linear stability of the achiral solution to chiral perturbations, and 
()64|) for the concentrations S2, /i2, -^^2 and M2, we obtain 



^5) 



Stability is determined by the eigenvalues of the matrix, which satisfy the cubic equation 

= + (2x + 5)A2 + ix' + 7x + 15/x)A + - |x/(9 + 30x). (86) 

Applying the Routh Hourwitz criteria (see the end of Section IXTl for details), we find that an 
instability occurs if / > 4^/(9 + 30%) ~ 2/15 (for large x)- 



4.2.3 Stability of achiral solution in Timescale II, smaller x 

At the end of Timescale II if x is small then the matrix in equation 



has the form 



d 

dt2 



e 



|a/i2 



I M2 M2 



2kfS2M2 (X-1)JV2 

M2 



a^2 



J 



M2 
X 

2kS2M2 



a^2 



e 



^7) 



All elements in the bottom row of this matrix are divergent since fi2, S2 ^ as t2 —>■ 00. 
When t2 becomes large the eigenvalues of this matrix are given by solutions of the cubic 







2kSoM 



2 X2 



2kSoMo 



N2_ axi^2N2 _ N2_' 
M2 2kS2M2 V M2, 



A 



2kfS2xN2 ^ x(x- l)Ar, 



M2^2 



(88) 



Applying the Routh-Hourwitz criteria for the signs of the real parts of the eigenvalues, we find 
stability of the achiral solution requires A > (which always holds), and C > which fails 
when / > aN 21^12/ '^kM2S2'i and a similar but more stringent inequality from the condition 
AB > C. So if 5*2 decays faster than ^2 then the symmetric solution is stable, and if 
S'2 and jj,2 decay at the same rate the stability depends on A^2 and M2 and requires / > 
l/i2{l-kM2/aN2)). 



4.2.4 Stability of achiral solution in Timescale III, smaller x 

With the scalings of ()76|) the matrix (jHHj) takes the form 
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N3 2/5 


A^s c-2/5 












2kfS3Ms 




2fcS3A/3 




M3 





\ 






e 


/ 





^9) 
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Note that the final row of entries are asymptotically larger than the other entries. Taking 
only the leading order entries for each term, the cubic governing stability can be written 



2kS^M^ W3 \ mJ) 2kS3Mi M3 ■ 



Rewriting this as A'^ + AX^ + BX + C, the Routh-Hourwitz criterion A > is always satisfied. 
The condition C > implies / < afi3N3/2kS3M3 which, at large times, reduces to / < 1, 
and as such is met in the large-time limit. Finally, AB > C fails if / > N^/lN^ + x^s) 
indicating an instability; whilst this depends on the unknown number-weighted and mass- 
weighted polymer concentrations, an instability is certainly possible for large enough values 
of the fidelity parameter, /. 



4.2.5 Stability of achiral solution in Timescale IV, smaller x 



Since the scalings for the concentrations S, fi, N and M are as in timescale III, the only 
change from (jSUI) is in the rate of growth of the perturbations since d/dt^ and d/dt^ scales 
differently with time and A^4, M4 are now time-dependent; we now have 
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dti 
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I Mi 





2kfS4MiE 
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N4, 
Mi 
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04^1 



e 



(91) 



As the solution progresses through timescale IV, the criterion for an instability to exist, 
namely / > N^/i^N^ + x^i)^ becomes easier to satisfy, since V4 decays slightly faster than 
M4, as shown by equations (jHT|) - (jH^ . 



4.3 Summary 

In this section we have analysed the kinetics of the concentrations of monomer, source and 
polymer as they evolve from zero to steady-state following a symmetric (achiral) solution. 

For larger values of x we have found three regimes through which the system evolves. 
Firstly there is a long induction period (of 0{e~^^^)) during which the source material builds 
up (until S becomes 0{€~^^^)); during this time the monomer concentrations remain small 
{0{e'^^^)) and the concentrations of polymer are extremely small {0{e)). This timescale ends 
abruptly as the catalytic feedback of polymer accelerates the breakdown of S into monomer. 
In the second timescale, which is very brief ((9(e^/^)), all monomer and polymer concentrations 
become large {0{e~^^^)) and decay towards steady-state. Finally, over the third timescale all 
concentrations converge to their steady-state values. 

In section 14.21 we analysed the linear stability of the growing achiral state through the 
sequence of timescales. By linear stability we mean that small external random forces which 
cause a chiral imbalance will be damped and reduce in amplitude. An instability indicates 
that such a perturbation will grow and so the system will undergo a symmetry-breaking 
bifurcation to occur during the evolution. 

As the system approaches steady-state in the third timescale we expect to regain the 
stability criteria / > fdx) with given by equation ()42p . However, this criterion for the 
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instability of the achiral solution is not valid in all the timescales; in the first timescale we 
find the alternative criteria of / > | for symmetry-breaking to occur. Note that this is 
independent of since in the first timescale the polymer and monomer concentrations are 
so low that inhibition of a homochiral polymer by a monomer of the opposite chirality is 
negligible; this leads to a considerable simplification of the linear stability analysis during the 
first timescale. In the second timescale we find an instability for / < 2/3. 

For / < 2/9 the achiral solution is always stable to such perturbations; whereas for 
2/9 < / < 1/2 the system is unstable to such perturbations only in the final stage of the 
kinetics, when the equilibrium solution is being approached and even then only for some 
values of x- For 1/2 < / < 2/3 the kinetics are more complicated since the system is unstable 
during the first timescale, but linearly stable during the second timescale and unstable during 
the final approach to equilibrium. Thus for these /-values the chiral purity may oscillate, but 
the system will eventually approach a chiral state. 

For smaller values of x fhe kinetics are more complex: there is still a long induction 
time followed by a period of rapid kinetics. This is more complicated, being split into vari- 
ous timescales; however the end result is always low concentrations of monomer and source 
species and large concentrations of polymer. There follows a long timescale over which the 
polymer concentrations reduce towards steady-state and the monomer concentrations and 
source species increase to steady-state. 

For smaller values of x ^i^n instability occurs in Timescale I for / > i, the instability 
persists in timescale II, now depending on / > 1/(2(1 — kM2/aN2)). In timescale III and IV 
an instability requires / > N/{N + x^), and at steady-state / > 4/5 — 34x/25. Thus, for 
smaller x there is still the possibility of a symmetry-breaking bifurcation occurring, though 
it requires a larger fidelity parameter. 

Although there advantages in x being large if one is seeking a symmetry-breaking bifur- 
cation, it should be noted that smaller x has other advantages, in that it allows polymers to 
form in larger concentrations, and these persist for longer times. 

5 Perfect fidelity 

In the extreme case / = 1 the feedback mechanism breaks down the precursor species (5) into 
chirally pure monomers with unit probability. Instead of a chirally pure homochiral steady 
state, there is a bifurcation to a state in which 6 asymptotically approaches ±1. We refer to 
such a state as a fully-bifurcated state. From equation (jH^ we see that such a state cannot 
arise if / < 1. The large-time asymptotics of the fully-bifurcated state differ significantly from 
/ < 1 since now there is no steady-state solution; instead we find unlimited growth of one set 
of homochiral polymer sequences and decay to zero for the sequences of opposite handedness. 
In this case (still ignoring e) the large time asymptotics are given by 

M = Mt, N = W/^, n = fit-^'^, S = St-\ (92) 

5 = 1-5, 6 = 1 ~e, r] = l-fj. (93) 

Assuming 5, 6 and r/ all decay to zero in the large time limit, we find the leading order 
equations 

So = kSM, M = a0, = a/l'^So = a0, (94) 
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S = -aN{x - lY^'^^O, 9 = f] = -axfif-^^^V. (95) 

These imply 

1 c<2/3 ol/3 

^ = p ^ = ^0, N=^, fl=^. (96) 

S ~ exp (-|aiV(x - l)t'/') , ^, ~ exp (-|ax^'/') • (97) 
Thus in terms of the concentrations of each chirahty we have 



i?l ~ /X ~ ^ , Rr^ N {^] , Pr^ M Sot, (9^ 



together with S" ~ l/H, and 

1^1 ~ L^t-^/3exp(-|(x-l)(a5ot2)2/3), (99) 
{L,Q} ~ {Z,g}t^/^exp(-|x(a5ot')'/'), (100) 

for some constants Li, L, Q. 

So in this case no finite steady-state solution is approached. As one might expect from 
the asymptotic expansions (03)^(03) i'^ ^^e case 5 — *■ ±1 we observe the unbounded growth of 
one type of homochiral sequence and, specifically, unbounded growth in the number of chains 
(iV), the mass of material in polymeric form (M) and the average length {M/N ~ (aSot'^Y^^). 
Concentrations of the sequences of opposite homochirality decay rapidly, and we expect that 
the average chain length approaches two, implying Q = 2L. 

All the above analysis has been for the simplified case for which e = 0; however, with one 
chain type decaying to arbitrarily small concentrations we may expect that the 0{e) term in 
is no longer negligible in this limit. Retaining the 0{e) term in this equation yields a 
slightly different scaling in the large time asymptotics for the high-fidelity case / = 1. We 
now have 

Ri-(AT, R-(^-^Y\ P-Sot, (101) 



.3aH 

together with Sq 1/kt and 

e £2 
Li ~ L ~ fl02) 

' kxiSa^Sif^y/^' k^X^So{3a^S^y/H^' ^ ' 

^ ^ Px'(3a252) 1/3^3- (103) 

Thus once again we see the less common homochiral polymer sequences assuming concentra- 
tions which decay to zero, albeit now with the simpler form of algebraic decay, the typical 
polymeric length again asymptoting to two (Q/L). The dominant homochiral sequences grow 
in number, mass and average length (with P/R ~ (a^ot^/S)^/^). As t — oo the chiral purity 
of the system approaches unity according to 

5 ~ 1 - r^, (104) 

kxSot 
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^" A;V52(9a5ot^°)^/=^' ^ ^ ~ ~^^^dW^^Wf^^' ^^^^^ 
Thus, we see the convergence to full chiral purity is more rapid for polymers than for 
monomers. 

In all the above analysis there is another chiral solution in which the left-handed ho- 
mochiral polymer sequences are dominant, and the right-handed homochiral sequences have 
concentrations which decay to zero asymptotically. 

6 Discussion 

After introducing our model in Section El we analysed its steady-states, and found that 
the symmetric steady-state solution exists for all parameter values, but that there are other 
solutions when the relative inhibition rate x and the fidelity / are large enough. The critical 
combination is 

(4 + 2x)(l + 2x) ,,„„. 
^ > ^= = (5 + 6x)(l+3x-)- 

At the point f = fc there is a supercritical pitchfork bifurcation, where two unstable steady- 
state solutions connect to the solution 6 = and make 5 = unstable for / > fc- When 
this inequality is satisfied there are two stable steady-state solutions with 6 non-zero, that is 
there are chiral solutions as well as the achiral solution, and the achiral solution is unstable, 
so any physical system will generically be attracted to one or other of the chiral solutions. An 
important effect to note from this formula is the role that cross-inhibition plays in making 
the asymmetric solutions accessible at low values of the fidelity. For small the fidelity has 
to exceed fc = 0.8 in order to obtain a symmetry-breaking solution; whereas at large this 
bifurcation point reduces to fc = 0.22 - a dramatic reduction. 

In Section^we analysed the kinetics of chain growth in a symmetric system, and found that 
there is a long induction time, during which a large stock of precursor chemical accumulates; 
an approximate, linear stability calculation shows that during this time, the achiral solution is 
unstable if / > 1/2. This behaviour is followed by a short timescale over which the precursor 
species is converted to monomers which are then polymerised. For this short time, monomers 
and chains are present in large concentrations. The concentrations of chains, monomers 
and precursor then all decay to their steady-state values, which, if the parameters /, x are 
such that an asymmetric steady-state exists, and the system has experienced some external 
perturbation away from the symmetric state, will be the chiral steady-state discussed earlier 
f section l3.2j) . In such a state both monomers and chains have a net chirality or handedness. 
Even for quite modest values of the chiral purity of monomer (say 6 = 0.7), the chiral purity 
of chains is extremely close to unity (6* = 0.990, t] = 0.995 at x = 2); see Figure 01 and 
compare Figures |3 and El 

Finally we have described the large-time asymptotics of the 'fully' bifurcated case which 
arises when / = 1, wherein chiral purities [6, 6 and rf) approach unity in the large time limit; 
this remains true even when we reintroduce the term which describes the slow spontaneous 
achiral decay of precursor species into both enantiomeric forms of monomer. 
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7 Conclusions 



In previous work we showed how quahtatively similar instabihties can lead to the massive 
amplification of self-replicating RNA polymer sequences over less efficient replicators (Wattis 
& Coveney 1999). 

In the present paper we have shown that an initially achiral system capable of stepwise 
polymerisation to homochiral polymer sequences with inhibition from the opposite-handed 
monomer, is subject to strong instabilities that drive the system overwhelmingly to one or 
other handedness for all homochiral sequences present. Mechanisms within this class may 
have played a role during the early stages of molecular evolution in determining the chirality 
of biologically relevant macromolecules, such as nucleotides and proteins, and there is ex- 
perimental evidence of this behaviour in the literature, for examples see Hitz et al. (2001, 
2002, 2003) and Joshi (2000). Although in the system studied by Joshi, addition of the cor- 
rect enantiomer to a growing polymer chain is more favourable than the wrong one, we have 
shown that if this cross-inhibition is stronger, then the system is more likely to undergo a 
symmetry-breaking bifurcation. 

Studies of this kind confirm the scope and power of modern methods of theoretical analysis 
for nonlinear dynamical systems of the kind that abound along the pathway towards the 
origins of life (Coveney, 1994). 
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